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Abstract 

We study a family of "classical" orthogonal polynomials which satisfy (apart from a 3-term recurrence 
relation) an eigenvalue problem with a differential operator of Dunkl-type. These polynomials can be ob- 
tained from the little q-Jacobi polynomials in the limit q = — 1. We also show that these polynomials provide 
a nontrivial realization of the Askey- Wilson algebra for q = — 1. 
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1. Introduction 



The Askey scheme [17], [18] provides a list of all known "classical" orthogonal polynomials. The term "classical" 
means that apart from a three-term recurrence relation 

P n +l(x) + b n P n (x) + UnPn-^x) = xP n ( X ), P_i = 0, P = 1 (1.1) 

(which is a general feature of orthogonal polynomials [8], [23]) these polynomials satisfy also an eigenvalue 
equation 

LP n (x) = X n P n (x) (1.2) 

where the operator L acts on the variable x. The operator L can be either a second-order differential operator 
(in this case we obtain the purely classical polynomials: Jacobi, Laguerre, Hermite) or a second-order difference 
operator on a uniform or non- uniform grid. The latter case leads to Wilson (q=l) and Askey- Wilson (g ^ 1) 
polynomials and their descendents. 

Bannai and Ito [4] proved a general theorem (which generalizes a famous theorem of Leonard [19], dealing 
only with polynomials orthogonal on a finite set of points) stating that all such orthogonal polynomials should 
coincide with the Askey- Wilson polynomials or their specializations. They also found a "missing" case of the 
Askey scheme corresponding to the limit q = — 1 of the q-Racah polynomials (see also [24]). The example of 
Bannai and Ito corresponds to polynomials orthogonal on a finite set of points. Hence in this case the operator 
L is merely a finite-dimensional 3-diagonal matrix which corresponds to the case of "Leonard pairs" (see, e.g. 
[24]). 

It is hence sensible to investigate other possibilities as q approaches -1 in the Askey scheme. Of course, the 
limit q = 1 is well studied and classified (see, e.g. [17]). The limit q = — 1 however, has not been explored much. 
In [1] Askey and Ismail have studied the limit q = — 1 for the g-ultraspherical polynomials, but in this case, the 
operator L disappears in the limit q = — 1. 

Here we show that there is a very simple class of polynomials which can be obtained from the little q- Jacobi 
polynomials in the limit q = — 1. Under appropriate choice of the parameters, the operator L survives in 
the limit q = — 1. The polynomials thus obtained are indeed classical: they satisfy the eigenvalue equation 
(1.2). But in contrast to the case of pure classical polynomials (like Jacobi polynomials), the operator L is 
a combination of a differential operator of first order and of the reflection operator R (see formula (2.14) in 
the next section). Operators of this type are known as Dunkl operators [11]. So far operators of Dunkl type 
were used typically to transform one polynomial family into another. In the one-dimensional case, the Dunkl 
operator is defined as [11], [6] 

T,f(x) = f'(x) + /(*)-/(-*) (L3) 

where \i is a deformation parameter (when fj, = the Dunkl operator becomes the ordinary derivative operator). 

Clearly, the Dunkl operator (1.3) reduces the degree of any polynomial by 1. Hence, there are no polynomials 
which are eigenfunctions of the operator T^. 
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Our operator (1.2), in contrast, preserves the linear space of polynomials of given degree. In the follow- 
ing we construct the general solution of the eigenvalue equation (2.15) and obtain explicit expression for the 
corresponding little -1 Jacobi polynomials P„ ^(x). 

We also show that the polynomials Pn~ X \x) are Dunkl-classical, i.e. T fJ/ Pn~ 1 \x) = Q n -i(x), where Q n {x) 
is another set of little -1 Jacobi polynomials (with different parameters) 



2. Limit of the little q- Jacobi polynomials as q — > — 1 

The little q- Jacobi polynomials are defined through the recurrence coefficients 

u n = A n —\C n , b n = A n + C n , (2-1) 

where A n , C n are given by 

_ „ (l-a g "+ 1 )(l- fl bg"+ 1 ) _ „ (1 - - bq n ) 

A n ~ 1 Tt 1 — 9 „ , i \ /-, , ,„,,, , — 0,q - 



(1 - abg2n+i)( 1 _ abq 2n + 2 )' * (1 - a6g 2 "+ 1 )(l - abq 2n ) ' 

They have the following simple expression in terms of the basic hypergeometric function 

Pn{x) = n n ( q n > ab( l n+1 q;qx ) (2 .2) 



aq 

with a normalization factor k„ to ensure that they are monic. 
They satisfy the orthogonality relation 

(bq; q)s 



J2 y^(aqyP n (q s )P m (q S ) = K S nm . (2.3) 
s =0 ( ' 9 ' 9 ' )s 



The moments corresponding to this weight function are 

(aq; q)n 



(2.4) 



(abq 2 ;q) n ' 

where (x; q) n = (1 — x)(l — qx) ... (1 — q n ~ 1 x) is standard notation for the g-shifted factorials (Pochhammer 
g-symbol) . 

There is a ^-difference equation of the form 

a(bq - x-^Pniqx) - P n {x)) + (1 - a;- 1 )(P„(g- 1 x) - P n {x)) = X n P n (x), (2.5) 

where 

\ n = (q- n -l)(l-abq n+1 ). (2.6) 

If a = q a , b = qP then in the limit q — > 1 we get the ordinary Jacobi polynomials with parameters a, (3. 
There is, however, another nontrivial limit if one puts 

q = -e £ , a = -e ea , b = -e e/3 (2.7) 
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and take the limit e — > 0. This is the limit q — — 1 of the little g-Jacobi polynomials. 
A direct calculation shows that in this limit, we have the recurrence coefficients 

_ (n+(l-0 n )a)(n + P + n a) _,_,y, (2n + l)g + ^ + a 2 + 

(2n + a + /3) 2 ' " ~ 1 > (2n + a + /?)(2n + 2 + a + /3) ' ^ 

where 

_ 1 + (-1)" 

On- g 

is the characteristic function of even numbers. 

The corresponding moments are obtained directly from the moments (2.4): 

(.iff ffiV "- 1 - 2 - 5 -- < 2 -»> 

where (x) n = x(x + 1) . . . {x + n — 1) is the ordinary Pochhammer symbol (shifted factorial). 

Using this explicit expression for the moments, we can recover the weight function w(x) for the resulting 
orthogonal polynomials. It is easily verified that 

w(x) = K \x\ a (l - x 2 )(' 3 - 1 )/ 2 (l + x), (2.10) 

where 

r(a/2 + 0/2 + 1) 

hi — 



r(/3/2 + l/2)r(a/2 + l/2)' 
Indeed, we have (using the ordinary Euler B-integral) 

J w{x)x n dx — c„, n = 0,1,2,..., 

where c„ is given by (2.9). The coefficient k is chosen to provide the standard normalization condition co = 1. 
Under the obvious conditions a > — 1, /3 > — 1, the weight function w(x) is positive, all moments c„ are well 
defined and the moment problem is positive definite, i.e. A„ > for all n = 0, 1, 2, . . .. 

Consider the form of the g-difference equation (2.5) in this limit. We divide both sides of (2.5) by e and 
introduce the operator L E which acts on any polynomial f(x) as 

L e f{x) = ae-\bq - x~ x ){f{qx) - f(xj) + ^(l - ar 1 )^" 1 ^ - f(xj) (2.11) 

(the parameters q, a, b depends on e as in (2.7)). For monomials / = x n we have in the limit e = 

L x n =i n x n +7 ln x n -\ (2.12) 

where 

C„ = 2(-l)" +1 n+(l-(-l)")( a + /3 + l), ^ = 2(-l) n n-(l-(-l) n )a. (2.13) 
This allows one to present the operator Lq in the form 

L = 2(1 - x)d x R + {a + (3 + 1 - aaT 1 )^ - ii), (2.14) 
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where R is the reflection operator Rf(x) = f(—x). 

Thus we have that our polynomials are classical: they satisfy the eigenvalue equation 



L P n (x) = X n P n (x) 



(2.15) 



where 




— 2n if n is even 
2(a + /3 + n+l) if n is odd. 



(2.16) 



But in contrast to true classical polynomial the operator L is not purely differential: it contains the reflection 
operator R. 

In the next section we construct the general solution of the eigenvalue problem (2.15) in terms of the Gauss 
hypergeometric functions. 

3. Different forms of ^-differential equation 

Consider the eigenvalue equation 



where Lq has the expression (2.14). The operator Lo is a differential operator of the first order containing 
the reflection operator R. Notice the important property of the operator L : it preserves any linear space 

of polynomials of degrees < N for any TV = 1,2,3, Hence the operator L behaves like the classical 

hypergeometric operator: for any n — 0,1,2,... there exists a polynomial eigenvalue solution L P n (x) = 
X n P n (x) for an appropriate set of eigenvalues X n . 

This solution can be found directly from the property (2.12) of the operator Lq. Indeed, in the basis x n the 
operator L is two-diagonal and hence the corresponding eigenvalue solutions can be found explicitly because 
the coefficients and rj n are known. 

We consider here another method, which, though more complicated, demonstrates the relation with the 
Gauss hypergeometric equation. 

We can arrive at a pure differential equation (without the operator R) by a standard procedure. Let us 
present the function F(x) as a superposition of the even and odd parts: 



where f(—x) = f(x), g{—x) = —g{x). The functions f(x) and g{x) are determined uniquely from F{x). The 
operator R acts on these functions as: Rf{x) — f(x), Rg{x) = —g(x). This allows one to rewrite the eigenvalue 
equation (3.1) in the form 



L F(x) = XF(x), 



(3.1) 



F(x)=f(x)+g(x), 



(1 - x)f'(x) + (x - l)g'(x) + (i + a + p- a/x)g(x) - X(f(x) + g(x))/2. 



(3.2) 



Consider also the associated eigenvalue equation 



RL F(x) = XRF(x) 



(3.3) 
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which is obtained from (3.1) by application of the operator R. This equation can be presented in the form 

(1 + x)f'(x) + (1 + x)g'{x) + (1 + a + + a/x)g(x) = X(g(x) - f(x))/2 (3.4) 
Adding and subtracting, we obtain the simpler system of equations 

f'(x) + xg'(x) + (l + a + P)g(x) = Xg(x)/2, 

xf(x) + g'(x) + ag(x)/x = -Xf(x)/2 (3.5) 

This casts the eigenvalue equation (3.2) equivalently as a system of two linear differential equations of first 
order. We can eliminate the function g(x) from this system: 

_ 2(x 2 -l)f(x)+Xxf(x) 
g[X) ~ 2(0 + 1) -A [A - b) 

and obtain a second-order differential equation for f(x): 

Ax{x 2 - l)f"{x) + 4((a + /3 + 3)x 2 - a)f{x) + \x(2(a + /3) + 4 - X)f(x) = 0. (3.7) 

By a change of variable x — > a; 2 one can reduce equation (3.7) to the Gauss hypergeometric equation thus 
obtaining the general solution: 

(A a+/3 i 1 _ A \ / A+2-2q 2/3+6-A \ 

4 ' 2 «±I ^J+^z 1 "^/ 4 3_ 4 ;x 2 j, (3.8) 

where Ci,C 2 are arbitrary constants. By construction, the function f(x) must be even for all values of the 
parameters a, p. This is possible only if C*2 = and the solution for f(x) contains only one undetermined 
constant: 

/ A a+i , 1 _ A \ 
f(x)=C 12 F 1 ( 4 ' 2 ^ 4 ;x 2 j. (3.9) 

From formula (3.6) we obtain the solution for the function g(x) 

A a+l _l 1 _ A 

2(a+l) 
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It is clear that the function g(x) is odd as desired. 

In formula (3.6) it was assumed that A ^ 2((3 + 1). If A = 2(/3 + 1), then the Gauss hypergeometric function 
2-Fi(x 2 ) reduces to the elementary function iF (x 2 ) and for the function f(x) we obtain 

f{x) = C l {l-x 2 )-^ L . (3.11) 

Similarly, if A = 2(/3 — 1) we obtain an elementary expression 



9{x) = ^ + 1 (1 - x 2 ) 2 



for the function g(x). 
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Now we have the general solution of equation (3.1) in the form 



A£M^!_a ^ AC(A) _ F (1 + |, + i - 
/ 2 (« + 1 ) 



F(x) = /(x)+ fl (x) = C(A) 2 F 1 4 ' 2 a ^ 4 ;* 2 --2^^ ' 4 ' 2 +3 ' 4 ;* 2 , (3.12) 



where the coefficient C(A) may depend on A (as well as the parameters a,0). 

We would like to get polynomial solutions, i.e. find eigenvalues A„, n — 0,1,2,... such that F n (x) is 
a polynomial in x of exact degree n. It is easily seen that polynomial solutions are possible only if either 
A = -An, n = 0, 1, 2, . . . or A = 2(a + (3 + 2 + 2n), n = 0, 1, 2, . . .. If A = -4n, the first term in (3.12) is a 
polynomial of degree 2n whereas the second term is a polynomial of degree 2n — 1. Hence for all A„ = — An 
we will have polynomials F n (x) of the even degree 2n. If A = 2(a + f3 + 2 + 2n) then the first term in 
(3.12) is a polynomial of degree 2n, while the second term is a polynomial of degree 2n + 1. Hence for all 
A„ = 2(a + (3 + 2 + 2n), n = 0, 1, 2, ... we obtain polynomials F n (x) having odd degree 2n + 1. This solves the 
problem, and we thus have the following explicit expressions. 

If n is even then 



If n is odd then 

P n 1 ( X ) = K n 



(_n n+a+P+2 \ / 1 _ « "+«+/3+2 



( l-n n+a+l+l \ I , o , , t\ / 1-n n+a+p+3 



(3.13) 



(3.14) 



where ft„ is an appropriate normalization factor to ensure that the polynomial Pn \x) = x n + 0(x n x ) is monic. 
(We need not its explicit expression). 

4. Relation with the symmetric Jacobi polynomials 

The symmetric Jacobi polynomials were introduced by Chihara [8]. They can be defined as follows. 
Let Pn ,r> ^ (x) be the (modified) Jacobi polynomials 

P^(x)^„ 2 F 1 (-^^|^ + 1 ;x) (4.1) 

which are orthogonal on the interval [0,1] 

Pi^ ) (x)P^(x)x^l - x)"dx = h n 5 nm 



I 

Jo 



>0 

(k„ is a factor needed for the polynomials P^' v (x) to be monic). 

We can introduce symmetric polynomials Sn'^ {x) by the formulas 



S&"Hx) = P^\x% S£l\(x) = xP£ +1 «\x 2 ). (4.2) 
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Then it is easily verified that the polynomials S n (x) satisfy the symmetric property S n (— x) = (—l) n S n (x) and 
are orthogonal on the interval [—1,1] 

J S^\x)S^(x)w(x)dx = h n S nm (4.3) 

with the weight function 

w [x) = \x\ 2i+1 {l-x 2 )^. (4.4) 

(Note that the original definition [8] by Chihara of the symmetric Jacobi polynomials S n (x) is slightly different. 
This is not essential for our purposes). 

If £ = —1/2, the polynomials sL 1 ^ 2 '' v \x) have the weight function w(x) = (1 — x 2 ) v and hence can be 
identified with the classical ultraspherical (Gegenbaucr) polynomials. This is why the polynomials S^ ,v (x) are 
sometimes called generalized Gegenbauer polynomials [5] . 

Starting from the polynomials S^ v (x) , we can perform the Christoffcl transform [22] 

» W - ^W; + ^»W , <4.5) 

where 

S&(-1) 
" ~ S*"(-l) ' 

The polynomials ^^(x) are again orthogonal on the interval [—1,1] with the weight function w(x) obtained 
from (4.4) as 

w(x) = w(x)(x + 1) = |x| 2 ^ +1 (l - x 2 )^{l + x) (4.6) 

(we do not take into account the normalization condition which is not essential for our purposes). 

The Christoffcl transform (4.5) can be presented in explicit form because the coefficients A n are expressible 
in terms of the hypergeometric function 2F\(z) with argument z = 1 (this is seen from formulas (4.2) and (4.1)). 

Hence from (4.5), we have an explicit expression for the polynomials S n (x) as a linear combination of two 
hypergeometric functions. The polynomials S n (x) are not symmetric, i.e. they satisfy the general 3-term 
recurrence relation 

S n+1 (x) + b n S n (x) + u n S n -i(x) = xS n (x), (4.7) 

where the coefficients b n , u n can be expressed in terms of the coefficients A n using the properties of the Christoffcl 
transform [22]. Hence these coefficients are also explicit. 

One can express the weight function (4.6) in an equivalent form 

w(x) = M 2 « +1 (l - x 2 )" +1 (l - x)- 1 

which corresponds to the Geronimus transformation [13], [26]: 

S n (x) = Si« +1 (x)-B n St«t\x) (4.8) 
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with some coefficients B n . 

We thus have that the same polynomials S n (x) can be obtained from the generalized Gegenbauer polynomials 
Sf l ' v (x) by either Christoffel or Geronimus transformations. 

In [9] these polynomials were presented in the form (4.8); the recurrence coefficients b n , u n were derived as 
well in [9]. 

The comparison of the weight functions (4.6) and (2.10) leads to the conclusion that the little -1 Jacobi 
polynomials P^ 1 \x) coincide with polynomials Sn' v \ where 

, a-1 0-1 

We thus identified our "classical" polynomials Pi~ X \x) with the " nonsymmetric" generalized Gegenbauer pro- 
posed by L.Chihara and T.Chihara [9]. 
Explicitly we have from (4.8) 

P^\x) = S^'^\x) - BnS^'^ix), (4.9) 

where 

_ 2rc+(l -(-!)*> 
2(a + /3 + 2n) 

The same polynomials were also considered in [2] , [3] from another point of view. 

We would like to stress that other basic properties of these polynomials , e.g. the existence of a Dunkl-type 
operator L providing the eigenvalue problem (2.15), as well as the origin of these polynomials as a limit case 
(<7=-l) of the little q- Jacobi polynomials had not been identified. 

5. The Dunkl-classical property 

All " classical" orthogonal polynomials satisfy an important characteristic condition: they are " covariant" with 
respect to a "derivative" operator V: 

VP n (x) = [n]P n -i(x), (5.1) 

where [n] is a specific function of n depending on the choice of the operator V, P n (x) is another set of "classical" 
orthogonal polynomials, and the operator D possesses the basic property of reducing the degree any polynomial 
by one. 

Well known examples of the operator V are: 

(i) the derivative operator V — d x ; 

(ii) the difference operator Vf(x) = f{x+ 1) — f(x); 

(iii) the g-derivative operator: Vf(x) = ^^zrj^r^ 

(iv) the Askey- Wilson operator Vf(x(s)) = £Mg^H£M±i^ , 

In the last case the function x(s) is either trigonometric x(s) = aiq s + a 2 q~ s + a or quadratic x(s) = 
a,2S 2 + a\s + a with some constants a*. 
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Recently it was recognized that apart from these operators there is one more operator which generates 
"classical" orthogonal polynomials. This operator is the Dunkl operator T u (1.3). Namely, in [6] it was shown 
that the only symmetric orthogonal polynomials P n (x) satisfying the property 

T„P n (x) = [n]„P n _i(x), [n]„ = n + (1 - (-1)> (5.2) 

are the generalized Hermite or the generalized Gegenbauer polynomials. Recall that symmetric orthogonal 
polynomials are defined by the property P n (— x) = (— l) n P n {x). The generalized Hermite polynomials H n ^ (x) 
[8] are symmetric orthogonal polynomials which are orthogonal on whole real line with the weight function 

w(x) = |x| 2M cxp(-a; 2 ) 

When fi = (i.e. in the case when the Dunkl operator T M becomes the derivative operator d x ) the generalized 
Hermite polynomials become the ordinary Hermite polynomials. 

The generalized Gegenbauer polynomials Sn' n \x) [8], [5] are orthogonal on the interval [—1,1] with the 
weight function (4.4). The generalized Gegenbauer polynomials satisfy the Dunkl-classical property (5.2) with 

M = £ + 1/2 

In both cases it is assumed that /i > —1/2. 

In the present case we have correspondingly the following simple but important result 

Proposition 1 The little -1 Jacobi polynomials Pn *\x) satisfy the Dunkl-classical property (5.2) with = 
a/2, where the polynomials P n {x) are again little -1 Jacobi polynomials with parameters (a, /3 + 2). 

The proof of this proposition follows easily from the explicit formula (4.9) and from the fact that the generalized 
Gegenbauer polynomials Sn' v \x) satisfy the Dunkl-classical propertry (5.2) [6]: 

T^\x) = [n],sH +1) (x), ^ = € + 1/2 

In contrast to the assumptions of [6], the little -1 Jacobi polynomials are not symmetric. Hence we perhaps ob- 
tained the first example of Dunkl-classical orthogonal polynomials beyond the family of symmetric polynomials. 
The problem of finding all such orthogonal polynomials is an interesting open question. 

All known families of "classical" orthogonal polynomials possess not only lowering operators like (5.1) but 
also raising operators & with the property 

<3P n (x) = u n+ iQ n+1 (x), (5.3) 

where the polynomials Q n (x) belong to the same family of classical orthogonal polynomials (albeit with different 
parameters). 

In the case of the little -1 Jacobi polynomials it is directly verified that the operator 6 does exists and has 
the expression 

Qf(x) = (x 2 - l)f'{x) + a{X ~ 1)2 f(-x) + (C8 + a/2)x - 1 - ^) f(x) (5.4) 
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Given (5.4), property (5.3) holds with 



1 - (-1)" 



and Q n (x) the same monic little -1 Jacobi polynomials with parameters (a, (3 — 2) . 

The generalized Hermite and Gegenbauer polynomials can be obtained from the ordinary Hermite and 
Gegenbauer polynomials through the acting of the Dunkl intertwining operator [11], [6]. 

Recall that the Dunkl intertwining operator acts on the space of polynomials by the formulas [11], [12] 

t t _ _ _ (1/2) n 

V^x =a n x , cr 2 „-i = <7 2n = 7 — -tttt; > 

(/i + l/2)„ 

and is realized by the following integral representation [11] 



i»r(i/2) 

It preserves the space of polynomials and has the fundamental intertwining property 

= V^d x . (5.5) 

From this property it is possible to obtain the following result (see Proposition 2 below). Assume that the 
monic polynomials P n {x) and Q n (x) are related as 

P' n {x) =nQ n - 1 (x). (5.6) 

Let us construct the monic polynomials P n (x) = cr^ 1 V^P n (x) 7 Q n (x) = o~ x V^Q n (x). Then these polynomials 
are correspondingly related: 

T^P n (x) = [n]„Q n _i(x). (5.7) 

In particular, if all polynomials P n (x),Q n (x), P n (x), Q n (x) are orthogonal then the operator allows to obtain 
Dunkl-classical polynomials (defined by property (5.7)) from ordinary classical polynomials (defined by property 
(5.6)). 

In [6] it was shown that the generalized Gegenbauer polynomials Sn' v ^ (x) can be obtained from the ordinary 
Gegenbauer polynomials S„ 1 ^ 2 ' v+fl \x) by the action of the intertwining operator V^: 

S^(x)=a- 1 V^- 1 / 2 ^(x), ,x = £ + l/2. (5.8) 

(A similar property for the generalized Hermite polynomials was obtained by Dunkl [11], [12]). 
Introduce now the ordinary monic Jacobi polynomials Pn^ (x) by the formula 

ptt*) (x) 2"(e + l) m (-n,n + ^ + V + ll_x\ 

Notice that this definition differs from (4.1) by an affine transformation of the argument. The polynomials (5.9) 
coincide with standard Jacobi polynomials orthogonal on the interval [—1, 1] [17]. 
We have the following 
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Proposition 2 The little -1 Jacobi polynomials (3.13), (3.14) can be obtained from the Jacobi polynomials 
(5.9) by the action of the Dunkl intertwining operator 

P^\x) = a~XPi U+1) (^ f = I* = f • (5-10) 

The proof of this proposition is based on formula (4.9) and property (5.8). 

6. Askey- Wilson algebra relations for the little -1 Jacobi polynomi- 
als 

The Askey- Wilson polynomials are closely related with the so-called AW(3)-algebra [25], [24]. Namely, there 
is a "canonical" representation of the AVF(3) algebra such that one of its generators has the Askey- Wilson 
polynomials as eigenfunctions. 

AVF(3) algebra usually consists of 3 generators [25]. Among different equivalent forms of the AVF(3) algebra 
we choose the following one, which possesses an obvious symmetry with respect to all 3 generators (see, e.g. 
[16]): 

XY - qYX = fi 3 Z + u 3 , YZ -qZY = mX + u u ZX - qXZ = fj, 2 Y + uj 2 (6.1) 

Here q is a fixed parameter corresponding to the "base" parameter in g-hypergeometric functions for the Askey- 
Wilson polynomials [17]. The pairs of operators (X, Y), (Y, Z) and (Z, X) play the role of "Leonard pairs" (see 
[24], [16]). 

The Casimir operator 

Q = (q 2 - l)XYZ + mX 2 + ll 2 q 2 Y 2 + fl 3 Z 2 + (q+ l)(wiX + uo 2 qY + UJ 3 Z) (6.2) 

commutes with all operators X, Y, Z. 

The constants u>j, i = 1,2,2 (together with the value of the Casimir operator Q) define representations of 
the ^W(3) algebra (see [25] for details). 

In the case of the little q- Jacobi operator, the realization of the AW(3) algebra is given by the operators 

X = g(L + l + qab), Y = x (6.3) 

where 

1 

9 ~ (q 2 - l)Va~b 

(the arithmetic meaning of the square root is assumed), and the operator L coincides with the difference 
eigenvalue operator for the little q- Jacobi polynomials in lhs of (2.5), i.e. 

Lf(x) = a(bq - x-^ifiqx) - f{x)) + (1 - x^^x) - f(x)) (6.4) 
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The operator Z is defined as 

Z/(a) = U-*)/(*/«) 



qy/ab 

We then have the AW (3) relations 

XY -qYX = Z + oj 3 , YZ-qZY = 0, ZX — qXZ = Y + uj 2 , (6.5) 

where 

6+1 o+l 

w 2 = -77 — TT; w 3 = ~" 



The Casimir operator 

g = (q 2 - l)XYZ + q 2 Y 2 + Z 2 + (q + l)(u 2 qY + u 3 Z) (6.6) 

takes the value 

Q = -b~\ 

These relations survive in the limit q = — 1. 
Indeed, let us define the operators 

X=^- 1 -±^±l, Y = x, Z=( X -1)R, (6.7) 

where Lq is the operator defined by (2.14) and the operator Y is multiplication by x. 
Then it is elementary to verify that the operators X, Y, Z satisfy the relations 

XY + YX = Z + a, YZ + ZY = 0, ZX + XZ = Y + f3 (6.8) 

which corresponds to the AW(3) algebra with parameters q = —1, w 3 = a, coi = 0, u; 2 = fl- 
it is easily verified that the Casimir operator commuting with X, Y, Z is 

Q = Y 2 + Z 2 . (6.9) 

In the case of of the realization (6.7) of the operators X, Y, Z, the Casimir operator becomes the identity 
operator: 

Q = I. (6.10) 

Note that relations (6.8) can be considered as an anticommutator version of some Lie algebra. Like in the 
case of an ordinary Lie algebra, the Casimir operator is quadratic in the operators Y, Z (the cubic part XYZ 
disappears in both "classical" limits g = ±1). 

The "canonical" representation of the algebra (6.8) is obtained in the basis P n (x) of orthogonal polynomials. 
In this basis the operator X is diagonal 

XP n (x) = K + l + a + P P n (x) 
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while the operator Y is 3-diagonal 

YP n (x) = P n+ i(x) + b n P n (x) + !i„P„-i(i). 

The monomial basis M n = x n provides another convenient representation, where the lower- and upper- trian- 
gular operators X, Y are 

XM n = £„M„ + r?„M n _i, YM n = M n+1 , 
with the coefficients £„, r\ n given by (2.13) . 

7. The special case a = and a "square root" of the Schrodinger 
operator 

Consider the special case a = of the little -1 Jacobi polynomials. From the expression for the weight function 
(2.10) it is seen that 

w(x) =k(1-x 2 )^- 1 ^ 2 (1+x). (7.1) 

The weight function (7.1) corresponds to the weight of the classical Jacobi polynomials [17] Pi a ' a+1 \x), where 
a = {(3 — l)/2. Moreover, from the Dunkl-classical property (5.2) it is seen that for a = the Dunkl operator 
T M becomes the ordinary derivative operator and we obtain the following realization of the Hahn property 

d x Pt a+1 \x) = nP£?' a+2 \x), (7.2) 

which characterizes classical orthogonal polynomials [15]. Thus, in the special case a — 0, the little -1 Jacobi 
polynomials become the ordinary Jacobi polynomials P^ L ' a+1 \x). For the sake of brevity, in what follows we 
will denote Q n (x) = P^ L ' a+1 \x) 1 i.e. Q n {x) are the little -1 Jacobi polynomials with a — and j3 = 2a + 1. 

It is well known, on the one hand, that the ordinary Jacobi P^ '^ (x) polynomials satisfy the Sturm-Liouville 
eigenvalue equation [17] 

SP^ b \x) = -\ n Pt h \x) 1 (7.3) 

where A„ = n(n + a + b + 1) and 

S = (1 - x 2 )d 2 x + (b-a-(a + b + 2)x) d x . (7.4) 

In the special case b = a + 1, we have the operator 

S = {l-x 2 )dl + {l-{2a + Z)x)d x . (7.5) 

On the other hand, we also know that the polynomials Q n (x) (being a special case of the little -1 Jacobi 
polynomials) satisfy the eigenvalue equation (2.15) 

L Q n {x) = X n Qn(x), (7.6) 
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where the operator L now is 

L = 2(l-x)d x R + 2(a + l)(l- R). (7.7) 

We thus obtain a rather surprising result: the ordinary Jacobi polynomials Pi a ' a+1 \x) satisfy not only the 
Sturm-Liouville equation (7.3) but also the Dunkl-type eigenvalue equation (7.6). 

Further analysis shows however that the operator L commutes with the operator So defined by (7.5): 

LoSo — SqLo = 

and that there is a simple algebraic relation between these operators: 

Ll - 4(1 + a)L = -45 . (7.8) 

Burchnall and Chaundy developed the theory of commuting ordinary differential operators V\ and Vi (the order 
of the operators may be arbitrary) [7]. They showed that if these operators commute, i.e. V{Vi = V 2 Vi then 
necessarily they satisfy an algebraic equation (see also [14]) 

V(V 1 ,V 2 )=0, (7.9) 

where V(V\, V2) is a polynomial of two such commuting operators. 

In our case we have a nontrivial example of the two commuting operators: So is the ordinary second-order 
Sturm-Liouville differential operator while L is a first-order differential operator of Dunkl type. These operator 
satisfy an algebraic equation (7.8) as well. 

There is also an interesting quantum mechanical interpretation of this example. 

Let us perform the change of independent variable x — siny in the operator S a : 

y \ cos^ y J 

Using the function 



= Vl + siny cos a+1 / 2 y, (7.10) 
we can transform this operator to the Schrodinger operator 

H=<£>(y)So<S>- 1 (y)=d 2 y -U(y) (7.11) 

with the potential 

U(y) = -(a + I) 2 + (a + 1/2)(a+ 2 1/2 - Slny) . (7.12) 

cos^ y 

This potential belongs to the class of generalized Poschl- Teller potentials and the corresponding Schrodinger 
equation can be solved explicitly (see, e.g. [10]). 

Under the change of variable x — siny and the similarity transformation L — > ^(y)L ^^ 1 (y) the operator 
L becomes 

Lo = 2 ( d y- ^oir) R y + 2 ( fl + x )' ( ? - 13 ) 
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where R y f(y) = f(—y) is the reflection operator with respect to the variable y. The constant term 2(a + 1) is 
inessential and we can introduce operator 



cosy 

without this term. Similarly, we can define the Hamiltonian 

Hl = -dj ]+ i* + VW* + V2-<**v) (7 . 15) 

y cos y 

without the superfluous constant term. (We changed the sign in the front of (7.11) in order to adjust to the 
canonical definition of a Schrodinger Hamiltonian [10]). 

It is then easily verified that the operator L\ commutes with the Hamiltonian H\: L\H\ — H\L\ and that 

L\ = H\. (7.16) 

Relation (7.16) shows that the operator L\ defined by (7.14) can be considered as a "square root" of the 
Schrodinger Hamiltonian (7.15). 

In what follows we will assume that a > 1/2. The potential U(y) of the Schrodinger Hamiltonian (7.15) is 
an infinitely deep well with walls at x = ±ir/2. Clearly, lim y ^__|- 7r y2 U(y) = oo. Hence the spectrum of the 
corresponding Hamiltonian is purely discrete and infinite. The corresponding wave functions tpn(y) are square 
intcgrable on the interval [— 7r/2,7r/2] and are solutions of the Schrodinger equation 

H!tp n (y) = E n ip n (y), n = 0,1,2,... 

with the energy 

E n = (a + n+l) 2 . (7.17) 

Explicitly the wave functions are 

1 — sin y 



i>n{y) = i>a{y) 2F1 

where 



-n, n + 2a + 2 
a + 1 



(7.18) 



Mv) = = Vl + siny cos a+1 / 2 y (7.19) 

is wave function of the ground state. Note that ^(±7r/2) = as required for eigenvalue solutions in the infinitely 
deep well. 

The operator L\ acts on the eigenfunctions as 

LiMv) = + n + l)My) (7.20) 

which corresponds to relation (7.16). 

This example can be generalized in the following way. 
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Consider the Schrodinger Hamiltonian 

J y 

with a potential U(y), and the Dunkl-type operator 



H = -dl + U(y) (7.21) 



L = (d + x (y))R (7.22) 

with some function x(u)- 

When do the operators H and L commute? It is easily verified that the commutativity condition HL = LH 
is equivalent to the conditions 

2 X '(y) = U(y) - U(-y) (7.23) 

and 

X"(y) = u'(-y) + x(y)(U(y) - U(-y)). (7.24) 

From condition (7.23) it follows that x'iu) 1S an °dd function. If we assume that xiv) has a finite derivative for 
y = then necessarily %(y) should be an even function, i.e. x(— y) — x{y)- 
Condition (7.24) can then be rewritten as 

X "(y) - U\-y) + x(y)(U(y) - U{-y)) = U'(-y) + 2 X (y) X \y) (7.25) 

from which we find 

X '(y) = -U{-y) + x 2 (y)-C (7.26) 
with some constant C. Changing sign in y we have similarly 

-x'(y) = -u(y) + x 2 (y)-c (7.27) 

Subtracting (7.26) and (7.27) we arrive at (7.23) which means that conditions (7.26) and (7.27) are compatible 
with (7.23) and (7.24). Adding (7.26) and (7.27), we have 

2x 2 (y) = U(y) + U(-y) + 2C. (7.28) 

There is a remarkable connection with the factorization method of the Schrodinger equation [10], [21]. Indeed, 
assume that for the potential V(y) and a constant C there is a the representation 

H + C = —dy + V(y) + C = A^A (7.29) 

where 

A = d v + X (y), Ai = -d y + X (y) (7.30) 

The Schrodinger Hamiltonian H is then said to be factorizable in terms of the operators A\ A. 
The potential V(y) is related with with the " superpotential" x{y) by obvious fromula 

x 2 (y)-x'(y) = v( y ) + c (7.31) 
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A refactorization performed by permuting the operators , A leads to the new Hamiltonian 

H + C = —dy + V(y) + C — AA\ (7.32) 
where the new potential V(y) is defined as 

x 2 (y) + x'(y) = v(y) + c. (7.33) 

Comparing formulas (7.31) and (7.33) with (7.26) and (7.27), we see that one can put 

V(y) = U(-y), V(y) = U(y), 

i.e. the refactorization procedure simply leads to a change of sign in the argument of the potential: U {—y) — > 
U(y). Such potentials were studied by Spiridonov and are special cases of more general self-similar potentials 
(see [21], [20] for details). 

We thus see that for any smooth even function x(u) the Dunkl-type operator L (7.22) commutes with the 
Hamiltonian H (7.21), where the potential U(y) is defined by (7.27) (the constant C can be arbitrary, of course). 
It is simply verified that in this case we have 

l 2 = (d y + x(y))R(d y + x(y))R = (d y + x(y))(-d y + xiv)) = -d 2 y + U(y) + C = H + C (7.34) 

i.e. the commuting operators L and H satisfy a simple algebraic relation of the second order. 

Thus one can say that the operator L is a " square root" of the Schrodinger operator for potentials U (y) . Of 
course, such operator exists only for potentials U(y) satisfying condition (7.27) with some even function x(y)- 

The physical meaning of the operator L is very simple. Indeed, let ip(x) be an eigenfunction of the Hamil- 
tonian (7.21) 

H^(y) = E^(y) (7.35) 

corresponding to the energy E. Obviously, 

LV(2/) - (±)VE^(y). (7.36) 

Acting on (7.36) with the operator R we get 

RL^{y) = A^(y) - [±)^E^{-y), (7.37) 

where the operator A^ is defined by (7.30). We see that the operator RL is the Darboux operator which 
transforms the eigenfunction ip(y) into the function ip(—y) which is the eigenfunction of the Hamiltonian 

H = -d 2 y + U(-y). 

Thus the operator RL "flips" the eigenfunctions corresponding to the potentials U(y) and U(—y). 
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